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ABSTRACT 


;ln  this  paper  the  tise  evolution  of  weak  shock  waves  propagating 
within  a  fluid-gas  Mixture  is  considered.  The  model  uses  continuum  classical 
theory  to  describe  the  shock  waves  and  allows  for  the  relative  notion  of  the 
gas  bubbles  and  liquid. \  It  is  shown  that  for  a  dissipative  system  the  waves 
are  governed  by  Burger's  Equation.  This  corresponds  to  the  far-field  solution 
for  waves  of  small  amplitude  and  long  wavelength.  The  nodal  also  provides  a 
non-linear  description  of  how  the  compression  wave  forms  and  approaches  a 
steady  state.  /(*?<*>  _y 


1 .  INTRODUCTION 


The  transient  behaviour  occur ing  when  a  weak  shock  wave  propagates 
into  an  inviscid  and  incompressible  liquid  containing  saall  gas  bubbles  is 
examined  in  the  long-time  regime.  A  semi-infinite  shock-tube  filled  with  a 
liquid-gas  mixture  and  fitted  at  one  end  with  a  moveable  piston  is 
considered.  When  the  piston  is  set  into  motion  at  constant  velocity  (u  ),  a 
wave  propagates  into  the  liquid-gas  mixture  and  is  modified  by  dissipative 
processes  and  by  non-linear  effects  associated  with  convection.  The 
dissipative  processes  involved  in  this  analysis  are  viscosity  and  thermal 
conductivity.  The  effects  produced  by  these  processes  upon  the  shock-wave 
propagation  can  be  complicated  since  both  waveform  and  velocity  are  modified 
under  their  influence.  Non-linear  effects  also  modify  the  wave  form  and 
become  important  when  of  comparable  magnitude  to  the  dissipative  processes. 
This  situation  will  arise  if  viscosity  and  thermal  conductivity  cause  velocity 
gradients  to  decrease.  Thus  in  the  far  field,  we  could  expect  the  exact 
solution  to  approximate  a  steady  compression  wave.  Since  viscous  forces  exert 
their  greatest  effect  in  the  region  of  the  shock  transition,  this  region  will 
be  discussed  in  detail. 


The  linear  theory  (van  Leeuwen  (1984))  starts  with  the  Navier-Stokes 
equations  for  a  fluid  gas  mixture  and  expands  the  dependent  variables  using  a 
perturbation  series.  ,  By  retaining  terms  to  order  0(e),  the  equations  of 
mass,  momentum  and  energy  can  be  solved  for  both  the  near  field  and  far  field 
solutions.  In  the  case  of  the  far  field,  the  solutions  are  given  in 
Appendix  1 .  These  solutions  oan  be  shewn  to  break  down,  when  t  becomes 
large,  and  this  feature  can  be  illustrated  as  follows.  The  width  of  the 
shock  transition  cone  (Thompson  (1972),  Gilbarg  and  Paolueei  (1953))  1st 
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for  the  case  where  the  velocity,  as  given  by  hi ,  is  of  order  0(/t).  Hence  as 
t  increases  the  size  of  the  shock  transition  zone  also  increases.  It  is  well 
known  however  that  S  is  of  order  0<e-')  for  steady  weak  shock  waves, 
where  e  is  an  appropriate  perturbation  parameter;  thus  when  t  becomes  large 
compared  to  e~  ,  the  linearized  solution  yields  an  excessively  broad  shock 
wave  and  so  fails  to  describe  the  phenomena  adequately. 


In  the  present  paper  a  uniformly  satisfactory  far-field  solution  is 
sought  for  a  shock  wave  propagating  into  a  liquid -gas  mixture.  Using 
"stretched”  coordinates  and  singular  perturbation  methods  a  uniformly  valid 
solution  for  t  ♦  •  can  be  derived.  Conditions  at  infinity  and  elsewhere  are 
assumed  uniform. 


An  analogous  type  of  piston  problem  has  been  studied  and  solved  by 
Moran  and  Shen  (1965)  (Referred  subsequently  to  as  MS),  for  a  viscous  heat- 
conducting  gas,  using  continuum  theory.  They  found  that  the  linearized 
Navier-Stokes  equations  are  valid  for  times  as  large  as  0(e~2)  mean  free  time 
after  the  piston  is  set  in  motion,  while  at  larger  times  the  solution  is 
governed  by  Burger's  equation.  Crespo  (1969)  has  considered  the  propagation 
of  infinitesmally  small  sound  waves  in  a  liquid  containing  gas  bubbles,  taking 
into  account  the  relative  motion  of  the  bubbles  and  liquid.  Also  the 
structures  of  shock  waves  in  liquids  containing  gas  bubbles  have  been  examined 
by  Campbell  and  Pitcher  (1958)  and  van  Wijngaarden  (1970).  Hie  methodologies 
of  MS  and  Crespo  are  closely  followed  in  this  paper.  A  continuum  description 
of  the  liquid-gas  mixture  is  given  in  section  2  together  with  simplifying 
assumptions.  The  method  of  analysis  employed  is  to  construct  the  exact  form 
of  the  Navier-Stokes  (NS)  equations  for  the  conservation  of  mass,  energy  and 
momentum  in  a  one-dimensional  viscous  unsteady  flow,  together  with  boundary 
conditions  appropriate  to  the  piston  problem,  the  dependent  variables 
appearing  in  the  NS  equations  are  then  expanded  in  the  usual  form: 


,.(0) 


+  e* 


(1) 


♦  ,V2' 


and  terms  of 
steady-state 
effects  come 
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order  0(e  )  are  retained.  The  NS  equations  to  order  0(e)  adopt  a 
fora  in  the  wave  frame,  while  to  order  0(e2)  the  time  evolution 
into  play. 


2.  BASIC  EQUATIONS  AND  MODEL 


The  liquid  host-medium  is  assumed  to  be  impregnated  with  gas  bubbles 
of  fixed  size  and  homogeneous  distribution.  The  liquid  is  also  assumed 
incompressible  and  inviscid.  (The  gas  obeys  the  perfect  gas  law.)  The  state 
of  the  liquid  can  be  described  in  terms  of  its  velocity  u,,  pressure  p^, 
density  p.  and  temperature  T,;  for  the  gas  the  corresponding  quantities  are 
u  ,  P  ,  pl  and  T .  The  volume  fraction  of  the  bubbles  in  an  aerated  liquid 
with  Initial  volume  VQ  is  denoted  by  0,  so  that  (1-0)  is  the  volume  fraction 
of  the  fluid.  The  liquid  is  taken  to  have  a  large  heat  capacity  and  that 
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variations  in  pressure  and  temperature  at  the  equilibrium  state  can  be 
neglected,  ie.. 


T 


i 


P 


i 


0 

0 


at  equilibrium 


The  above  condition  with  respect  to  pressure  implies  that  the  radius  of  the 
gas  bubbles  is  small  compered  to  the  width  of  the  shock  transition  rone. 
Further  assumptions  made  are  (i)  that  the  viscous  forces  are  attributable  to 
gas  bubble  drag,  (ii)  that  heat  transfer  occurs  only  between  the  gas  bubbles 
and  surrounding  liquid,  and  (iii)  that  the  thermal  notion  of  the  gas  bubbles 
is  negligible  and  so  does  not  contribute  to  the  pressure. 


Both  dependent  and  independent  variables  are  ini£i£l\y^eypressed  in 
a  coordinate  system  designated  in  terms  of  aA*  (ie.,  f  •  (u,P,p,8,T)),  but  for 
convenience  the  *  is  dropped  from  all  equations. 


In  the  nomenclature  described  above ,  the  exact  Navier-Stokes 
equations  (see  Ockendon  and  Tayler  (1983))  expressing  the  conservation  of 
mass,  momentum  and  energy  of  a  liquid-gas  system  became  for  a  one-dimensional 
unsteady  flow  (van  Leeuwen  (1984),  Crespo  (1969),  MS)  as: 


(i)  Equation  of  continuity: 

3t(1-6)  +  3x(1-B)u4  -  0, 

(ii)  Equation  of  continuity  for  the  gas: 

3t{Ppg)  +  “  °* 

(iii)  Equation  of  momentum  for  the  liquid-gas  composition: 
Pt(1-8> (®tut  +  +  3xPi  “ 

(iv)  Equation  of  momentum  for  the  gas: 


T  9  4WP«  2/3 

3xPi  '  2  Wig  +  V*V  '  2  ptVi(-ir)  uig  “  °' 


Share  u. _  •  u.  -  u  . 

ig  4  g 


(2.1a) 


(2.1b) 


(2.1c) 


(2. Id) 
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(v)  Equation  of  energy  for  the  gas: 


3.P.  +  u  3  P4  - 
t  t  g  x  i  p 


9  * 


Op  +  u8p)  - 
t  g  g  x  g 


4*p  2/3 

3(lfg.  1)oA(^i)  (Tt  -  T  )  -  0  , 


(2.1e) 


where  3  denotes  partial  differentiation  (elsewhere  the  alternative  notation 
' will  also  be  used). 


(vi)  The  equation  of  state  for  a  thermally -perfect  single-component, 
dilute  gas  is  given  by 

pg  “  W  <2.1f) 


where  R  is  the  universal  gas  constant.  The  corresponding  liquid  equation  of 
state  is  simply  p^  »  constant. 

In  the  above  equations  v,  and  a.  are  the  kinematic  viscosity,  and 
the  thermal  conductivity  of  the  fluid  respectively,  v  is  the  ratio  of 
specific  heats  of  the  gas,  and  T  is  a  function  which  depends  on  ft  and  which 
for  6  small  is  essentially  unity  (Crespo  (1 969)). 


One  approach  to  solving  the  above  equations  is  to  apply  the 
fundamental  principle  of  dimensional  analysis  (ie.,  every  problem  must  be 
expressible  in  terms  of  dimensionless  variables).  Firstly  we  seek 
dimensionless  variables  for  the  independent  variables  (x,t).  An  examination 
of  the  physical  constants  appearing  in  (2.1)  reveals  shows  that  a  constant 
with  dimensions  of  length,  can  be  formed  which  renders  x  and  t 
dimensionless.  That  is,  we  can  choose  without  loss  of  generality  a  length 
scale  k  «  aereup1 •  where  re  is  the  radius  of  a  gas  bubble  and  a#  «  /rt^  is  the 
equilibrium  speed  of  sound  in  the  liquid,  such  that: 


x  a  kx,  (2.2a) 

t  -  fc@t.  (2.2b) 


Similarly  the  dependent  variables  9  *  (u,P,p,£,T)  can  also  be  reduced  to 
dimensionless  form,  following  the  example  of  MS. 


The  dimensionless  fluid  variables  can  be  expressed  asymptotically  in 
terms  of  fluid  variables  as  a  power  series  in  c  about  the  equilibrium 
state  {0,PB,P#,8B,Ta}  ie., 

♦  -  te<1  +  «t(1)  ♦  «V2)  ♦  •••),  (2.3a) 
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for  the  state  variables  P,p,$  and  T,  and  as 

♦  -  eae  +  e^*2*  +  ...),  (2.3b) 

A 

for  u.  The  dimensionless  perturbation  variable  e  is  chosen  as  a  measure  of 
the  piston  Mach  number 

u 

e  «  (2.4) 

e 

where  up  is  the  initial  piston  velocity. 


The  linear  theory  fails  When  the  diffusion  thickness  and  the 
compression  wave  thickness  are  equal,  ie.,  when  t  =  e-2.  To  find  a  uniformly 
valid  solution  for  t>e~2  (ie.,  for  the  far  field)  where  the  linear  theory  is 
no  longer  valid,  it  is  necessary  to  introduce  "stretched"  coordinates  for  the 
dimensionless  independent  variables,  in  order  to  take  into  account  slow 
variations  in  the  wave  form.  This  can  be  achieved  by  introducing  the 
following  scale  transformation  for  the  independent  variables  (Su  and  Gardiner 
(1969)): 


5  -  en(x  -  Xt),  (2.5a) 

t  -  en+1t,  (2.5b) 


where  X  is  a  dimensionless  constant  related  to  the  velocity  of  propagation  of 
the  disturbance,  and  is  to  be  determined  by  the  analysis.  The  quantity  n  is  a 
positive  constant  (also  to  be  determined)  such  that  time  variations  of  any 
state  variable  are  balanced  by  both  non-linear  and  dissipative  effects. 


In  the  present  nomenclature,  the  system  of  equations  (2.1)  can  be 
expressed  in  dimensionless  variables.  The  scheaw  is  illustrated  by 
considering  the  continuity  equation  (2.1a).  After  a  change  of  independent 
variables  from  (x,t)  to  (x,t)  equation  (2.1a)  becomes: 


ae3t(l-6)  +  3x(1-8)ujl  -  0. 


(2.6) 


employing  (2.5a)  and  (2.3)  we  obtain  to  order  0(e)  the  steady  state  equation: 

(2.7a) 


xeeV(1>  *  <1"8e)8|ul1)  "  °' 


and  to  order  0(c): 


*®eV<2)  '  Vt0<1)  *  (1“Pe)3Cut2)  “  V<1,V»1) 
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(2.10a> 


-  vi'V 


(1) 


o. 


which  new  incorporates  tine-evolution  effects.  Proceeding  in  a  like  aanner 
for  the  reaaining  equations  we  obtain  to  order  0(e): 


T(1)  =  u(1>  -uj1*  =  p(1>  -  p(1>  =  0, 

g  g  i  g  g 


3eug1)  -  *Vpg1}  +  9<1))  “  °* 


(2.7b) 

(2.7c) 

(2.7d) 


Integrating  equation  (2.7)  and  equating  the  constants  of  integration  to  zero 
gives  the  state  variables  of  the  gas  in  terns  of  the  fluid  velocity: 


ui”(O.X(1-8,). 


<1-@J 


(2.8) 


where  the  dimensionless  velocity  X  satisfies  the  relation 


B  (1-B)' 
e  e 


(2.9) 


together  with  the  obvious  physical  constraint  that  the  fluid-gas  mixture  is 
not  purely  fluid  (0e  -  0)  or  gaseous  (6@  *  1).  This  velocity  has  a  minimum 
value  when  0^-  1/2. 


^Following  a  similar  procedure  to  the  above  for  the  equation  of 
order  0(e/),  in  which  the  effects  of  dissipation  are  embedded  in  the  flow,  we 
find  that  the  parameter  n  should  be  chosen  as  unity.  Tor  this  choice  of  n  the 
second  order  perturbation  equation,  together  with  (2.7),  can  be  expressed  as: 


-  U ><2)  ♦  0(2)+  P(,)0(1>)  ♦  3 ,(u<2)  ♦  p^V1*  ♦ 

V  g  g  J  V  g  g  t  i  J 


+  3t(p^1)  +  eM))  -  o, 


(2.10b) 


X0e$l1)3euJ1)  -  X(1-ee>aeu^2)  +  (1-^)3^’ *  +  (1-e#)u{1)3eu[1) 


+  -  0, 


(2.10c) 
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(2.10d) 


Vg  +  ♦«“*  '  g  '  " 


X(1-Y  )3  P<n  -  T<2)  -  0, 

9  E  9  «pTT  g 


-  42'  -  t<21  -  0. 


(2.10e) 


(2.10f) 


The  constant  tt  »  4„0Pr  '  i®  ****  themal  relaxation  coefficient. 

f  e 


3.  BURGER'S  EQUATION  TOR  THE  SHOCK  WAVE  STRUCTURE 


First-order  perturbations  in  the  state  variables  appearing  in  (2.10) 
can  be  written  in  terns  of  5 ,  via  the  steady  state  equations  (2.7). 
Equations  (2.10)  can  then  be  solved  to  give.*'?'  in  terns  of 
Uj  |  and  Uj  .  The  partial  derivatives  f  ^  'are 


v<2> 


(3.1a) 


Vg2>  "  *2  + 


(3.1b) 


-  it  (».S  *  Vi2')’ 


(2)' 


(3.1c) 


3,u*2)  *  3Eui2)  -  1  (1"8e>3EEul1>' 
Eg  E  *  f  *  EE  * 


(3. Id) 


T<2)--^f  (1  .T  )u <1) 

9  Ve8e  9  4,5 


(3.1e) 


«here  G.  (j»i,2,3)  are  functions  of  the  first  order  perturbation  variables 


(3.2a) 


G1  =  3Te(1)  +  3c(u^1)6(1)), 

G2  =  (1-ee)3TuJ1)  +  3ir(uJ[1)<7(1-Be>'41>  +  Xee^'  (3.2b) 

G3  =lVi  Vg1’  +  Pg1)<Ui1)  ‘  ^(1>))*  (3*2c) 


Dif ferentating  equation  (3.2e)  with  respect  to  £,  and  substituting 
Tg2^  using  (2.10f)  together  with  (3.2),  we  find  that  u^1 '  satisfiei 


for 

satisfies  Burger's 


equation  (Whitham  (1974)): 


0  3  uj  1  * 

ext 


(1 ) 


2  3?CUt 


(1) 


=  0  for  t  >  0, 


(3.3) 


where 
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uPTT(ira~1 


is  a  dimensionless  constant  which  behaves  as  a 


viscosity-like  term  and  is  typically  the  inverse  of  the  Reynolds  number. 
This  term  serves  to  prevent  singularities  which  would  occur  if  9=0.  in 
equation  (3.3)  u^u^  ^  is  a  non-linear  convective  term  and  V2  9u^  is  a 
diffusive  equation  viscous  term. 


The  algebraic  form  of  the  above  equation  confirms  that  the  stretched 
coordinates  (2.5)  are  appropriate,  since  equation  (3.3)  retains  terms 
describing  time-dependent  flows  with  non-linear  convection  and  dissipation. 


4.  SOLUTIONS  FOR  A  DISCONTINUOUS  SHOCK 


The  general  initial-value  problem  for  Burger's  equation  (3.3)  can  be 
solved  exactly  by  means  of  a  Cole-Hopf  transformation  (Cole  (1951),  Hopf 
(1950)).  This  transformation  reduces  Burger's  equation  to  a  classical  heat- 
conduction  equation  which  can  then  be  solved  using  integral- transform 
methods.  In  a  coordinate  system  (I,t),  where  \  =  £  and  t  =  0@  t,  we  can 
introduce  a  function  £(£,t)»  which  is  related  to  u  ,  in  such  a  way  that: 


r  1  JUr 

"  "  9  U*  Z’ 


_  1  f1  (1)  6  (1K_ 

£,—  =  "T  [t  u.  —  —  U-  r  )  £  , 
T  6  v2  l  2  1,5/ 


(4.1a) 


(4.1b) 


and  this  maybe  achieved  by  noting  that  the  integrability  condition  is 
precisely  equation  (3.3)  (i.e.,  3-  (4.1a)-3=(4.1b)  ■  0,  yields  equation 
(3.3)).  This  non-linear  transformation  eliminates  the  non-linear  terms 
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n^-,- 


appearing  in  (3.3).  Integration  of  (4.1a)  then  yields: 

I  »  exp(-  —  Ju^dt).  (4.2) 

On  substituting  equation  (4.2)  into  equation  (3.3)  we  find  that  £  is  a 
solution  of  the  classical  heat  conduction  equation: 


The  initial-value  problem  for  Burger's  equation  is  therefore  exactly  soluble 
and  pseudo-stationary  solutions  can  be  obtained  from  the  ansatz: 


*(?,t)  -  *(£  -  at). 


from  which  it  follows  that 


*(£)  -  a(l  +  exp(-c( £  -  CQn). 


More  explicitly,  solutions  of  (4.3)  corresponding  to  a  discontinuous  shock 
wave  at  (  •  o  can  readily  be  found,  for  the  case  that  the  initial  condition 
at  t  =  0  is 


(5,0) 


uQ,  5  <  ° 
0  ,  I  >  0, 


equation  (3.2)  yields: 


£(£,0) 


uf 

exp(-  -g-) 


5  <  0 
£  >  0, 


where  £  is  a  continuous  function  through  £  ■  0  and  uQ  is  an  arbitrary  constant 
which  can  be  set  equal  to  unity  without  loss  of  generality.  This  initial 
condition  is  the  one  which  matches  the  asymptotic  form  of  the  near-field 
solution  (see  appendix  1).  The  solution  of  (3.3)  with  this  initial  condition 
can  be  found  by  taking  a  Laplace  transform  of  (3.3)  with  respect  to  t.  This 
yields,  in  untransformed  coordinates  (Ockendon  and  Tayler  (1983)): 


(lot) 


1-1 


1  + 


l-VLT 

/  29t 


(4.5a) 


» 


where  erfc(x)  *  l-erf(x)  “  1  “  hi  fo  e_<*  da*  The  initial  condition  can  be 
recovered  from  (4.5a)  if  £  t  0  and  t  =  0.  In  coordinates  (x,t)  the  fluid 
velocity  (4.5a)  nay  be  expressed  as 


u*1J (x,t) 
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(4.5b) 


This  solution  is  similar  to  that  of  MS  for  shock  waves  propagating  into  a 
viscous  heat-conducting  gas  in  the  far  field.  The  MS  solution  is  recovered 
if  0  ♦  28(y  +  1),  Be  -*■  2(y  +  1)_1  and  X  ♦  /y. 

If  we  let  t  ♦  «®  with  5  kept  finite,  then  (4.5b)  reduces  to 
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(4.6) 


where  the  following  properties  of  error  functions,  erfc(0)  =  1  and 
erfc(-  •)  =  2  have  been  used.  This  solution  represents  a  shock  wave  of 

u  e 

constant  profile  propagating  with  a  velocity  of  X  +  -3*—.  This  velocity  is  a 

pe 

second-order  approximation  to  the  Rankine-Hugoniot  result  for  the  velocity  of 
a  weak  shock  (cf  Taylor  (1910)).  the  first-order  approximation  is  given  by 
(2.9).  He  further  observe  that  if  the  solution  (4.6)  is  stressed  as 

g(x,t)  -  -  *  u*1J(x,t) 


then  g(x,t)  satisfies  the  steady-state  Burger  Equation  (cf  MS  where  the 
asynptotic  solution  is  nissing  a  factor  of  V, )  ie. , 


5.  NUMERICAL  ILLUSTRATION 


Hie  forgoing  solution  (4 .5b)  is  illustrated  for  a  shock  wave 
propagating  into  water  (p,  ■  103  kgm“3,  ae  *  1.51  x  103  as-1, 

■  8.57  x  10"  m2s-1 ,  a*  *6.1  x  10_1  J/(ms°X),  -  300*X)  containing  gas 

bubbles  of  radius  r„  -  10*  a.  In  Figure  1  and  4  we  consider  cases  where 
e  =  .05  ot  e  «  .2,  that  is  the  voluae  fraction  of  gas  in  the  fluid  is  either 
5%  or  20%.  The  velocity  profiles  are  presented  for  both  the  linear  (Appendix 
A-1)  and  non-linear  cases  in  diaensionless  units. 


For  a  small  time  interval  the  velocity  gradients  are  very  steep  and 
linear  theory  is  then  generally  in  good  agreement  with  the  non-linear 
theory.  This  result  is  expected  since,  for  t  small,  dissipative  terms  are 
more  important  than  non-linear  convective  terns.  As  t  increases,  the  wave 
propagates  further  into  the  liquid-gas  mixture  and  the  velocity  gradients 
diminish,  leading  to  a  broadening  of  the  shock-wave  i.e. ,  non-linear  effects 
begin  to  dominate. 


When  e  ♦  0  the  linear  and  non-linear  solutions  coalesce  and  shock 
fronts  become  infinitely  steep;  see  Figures  1  and  2.  Conversely, 
if  e  increases  then  the  two  solutions  diverge  and  the  velocity  gradients 
diminish. 


The  effects  of  altering  the  viscosity-like  term  8  are  shown  in 
Figures  5  and  6.  Increasing  8  (Figure  6)  leads  to  a  smearing  of  the  velocity 
profile,  which  might  be  expected  since  8  contains  terms  related  to  momentum 
and  thermal  relaxation.  If  8  is  reduced,  on  the  other  hand,  then  the  velocity 
gradient  becomes  very  large  (Figure  5). 


6.  DISCUSSIOM 


The  propagation  of  pressure  waves  into  gas  bubbles  suspended  within 
a  liquid  is  of  some  considerable  interest,  in  this  paper  a  mathematical  model 
is  developed  which  investigates  the  evolution  of  a  shock  wave  propagating 
within  a  two-phase  fluid,  and  proceeds  with  a  description  of  how  the 
compression  wave  forms  and  approaches  the  final  steady  state.  The  model 
developed  is  an  extension  of  the  previous  linear  model  of  van  Leeuwen  (1985) 
and  is  also  a  generalisation  of  the  work  of  Moran  and  Shen  (1966),  who 
considered  the  case  of  a  shock  wave  propagating  within  a  viscous  heat- 
conducting  qas.  The  linear  and  non-linear  solutions  of  the  flow  that  have 
been  developed  give  a  uniformly  satisfactory  description  of  the  flow. 

Important  features  of  the  present  model  are  that  the  Mavier-Stokes 
Equations  for  a  liquid-gas  mixture  are  solved  using  a  perturbation  expansion 
on  the  dependent  variables  and  a  dilation  of  the  independent  variables.  This 
procedure  extends  the  solution  beyond  the  region  where  the  linear  solution 
breaks  down.  Decoupling  the  resultant  partial  differential  equation  leads  to 
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Burger's  equation  for  the  velocity  of  the  liquid.  By  using  Laplace  transform 
techniques  this  equation  can  be  solved  in  closed  form. 


For  small  values  of  time  the  non-linear  solution  agrees  well  with 
the  linear  one,  however  for  the  far  field  the  liquid  velocity  in  the  non¬ 
linear  case  runs  ahead  of  the  linear  one. 
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APPENDIX  1 


the  solution  of  the  Navier-Stokes  equations  for  unsteady  flow  to 
order  0(e)  has  been  derived  for  a  shock  wave  propagating  into  a  liquid-gas 
■ixture  (van  Leeuwen,  1984) .  the  dimensionless  solution  for  large  tiaes 
(i.e.,  the  high  frequency  case)  concomitant  to  (3.5b)  can  be  written  as: 
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and  for  the  density,  pressure  and  teaperature  of  the  gas  as: 
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the  velocity  of  the  gas  is  simply: 


FIGURE  5 
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